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A necessary condition for superconductivity (SC) driven by electron correlations 
is that electron-electron (e-e) interactions enhance superconducting pair-pair corre¬ 
lations, relative to the noninteracting limit. We report high-precision numerical cal¬ 
culations of the ground state within the frustrated two-dimensional (2D) Hubbard 
Hamiltonian for a wide range of carrier concentration p (0 <P< 1) per site. We find 
that long range superconducting pair correlations are enhanced only for p ~ 0.5. At all 
other fillings e-e interactions mostly suppress pair correlations. The enhancement of 
pair correlations is driven by the strong tendency to local sin^et bond formation and 
spin gap (SG) in p = 0.5, in lattices with quantum fluctuation®^. We also report de- 
terminantal quantum Monte Carlo (DQMC) calculations that are in strong agreement 
with our ground state results. Our work provides a key ingredient to the mechanism of 
SC in the 2D organic charge-transfer solids (CTS), and many other unconventional su¬ 
perconductors with frustrated crystal lattices and p ~ 0.5, while explaining the absence 
of SC in structurally related materials with substantially different p. 


The possibility that e-e interactions can be the driving force behind SC in correlated-electron systems has been 
intensely investigated since the discovery of SC in the high T^ cuprates. The minimal requirements for a complete 
theory are, (i) the superconducting pair correlations are enhanced by e-e interactions, and (ii) pair correlations are 
long range. For moderate to large e-e interactions, pair correlations are perhaps best calculated numerically, which 
however can be done only for finite clusters. The simplest model incorporating e-e interactions is usually assumed to 
be the Hubbard model, which in quite general form can be written as 

H = — 'y ^ H.c.) -f U y ^ + — y ^ VijUiTij. (1) 

i (ij) 

All terms in Eq. [^have their standard meaning. The first sum is the kinetic energy of noninteracting electrons within 
a 2D tight-binding model with hopping matrix elements tij', U and Vij are onsite and nearest neighbor (n.n.) Coulomb 
interactions respectively. Existing numerical calculations within Eq. [l] on 2D lattices have failed to find enhancement 
of pair-pair correlations relative to the noninteracting model without making severe approximation^. 

It has sometimes been surmised that correlated-electron SC might evolve upon doping a spin-g app ed semiconductor, 
as would occur in toy models such as a 2D lattice consisting of weakly coupled even-leg laddera^EI. Finding realistic 
2D models with SG and enhanced pair correlations however remains challenging. In the present work we demonstrate 
from explicit numerical calculations on frustrated 2D lattices enhanced pair correlations evolving from a spin-gapped 
state at a carrier density p ~ 0.5, far from the region most heavily investigated until now (0.7 < p < 1.0). We 
further point out the strong relevance of the resulting theoretical picture to real materials, in particular the 2D CTS 
superconductors, which were discovered earlier than the cuprate^ but are still not understood. 

There occurs an effective e-e attraction uniquely at p = 0.5, in lattices with strong quantum fluctuations, driven by 
charge-spin-lattice coupling. Consider the four-atom dimerized “molecule” of Fig. [^a), with two strong intradimer 
bonds and one electron on each dimer. In the absence of the interdimer bond, electron populations are equal on all 
sites. As the interdimer electron hopping is switched on slowly, there is net migration of charge to the two center 
atoms, due to the attractive antiferromagnetic spin-coupling whi ch l eads to a singlet boncPI. The charge migration 
is enhanced by the presence of electron-phonon (e-p) interactional^. The effective attraction is stronger than that 
near p ~ 1, where the tendency to such charge migration is necessarily smaller, with the neighboring sites already 
occupied. The charge-ordering (CO) of Fig. ma) in the spin-singlet state persists in the thermodynamic limit in one 
dimension (ID) p = 0.5, where for V < Vc(Uj e-e and e-p couplings act co-operativelj^ to give the spin-Peierls state 
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with the 2ki;’ bond modulations of Fig. and (c). Phase segregation is avoided, as the key requirement for charge 
migration is a final state with vacancies on both sides of the singlet bond (see Fig. [^a)). The spin-Peierls state at 
p = 0.5 is a paired-electron crystal (PEC), in which singlet-coupled n.n. singly occupied sites are separated by pairs 
of vacancies. Similar PECs occur in the p = 0.5 zigzag ladder (Fig. and in the 2D anisotropic triangular lattice 

(Fig-B e)) for sufficiently large lattice frustratiorP. We have not found the PECP at any other p. The exceptional 
stability of the PEC at p = 0.5 is due to its commensurate structure. The PEC has been experimentally observed in 
a number of CTpl (see Supplemental Information). 

Based on a valence bond (VB) perspective that has some overlap with Anderson’s resonating valence boncPSl 
approach to the nearly p = 1 limit, we posit that SC is achieved in p ~ 0.5 upon destabilization of the PEC, either 
due to increased frustration or very weak doping. The PEC wavefunction is dominated by covalent VB diagrams 
with periodic arrangement of the n.n. singlet bonds. Close to the PEC, the static CO and bond order are lost, 
but we anticipate the wavefunction to continue to be dominated by VB diagrams with singlet bonds between n.n. 
charge-rich sites, except that the arrangement of the bonds is no longer periodic. One such diagram is shown in 
Fig. Bf)(i)- Within Eq. pairs of VB diagrams with only n.n. bonds are coupled through the diagrams with next 
nearest neighbor (n.n.n.) bonds, as in Figs. [^f)(ii) and (iii). We collectively refer to diagrams with only n.n. and 
n.n.n. bonds as those with short bonds. There will be considerable pair tunneling in a wavefunction dominated by 
VB diagrams with short bonds, and we will refer to such a wavefunction as a paired electron liquid (PEL). The PEL 
would be most stable near p = 0.5, where the gain in kinetic energy due to pair tunneling is the largest. 

We consider an anisotropic triangular lattice with tij = {t^, ty, tx+y}- Vij can similarly have three components. 
We express all quantities with dimensions of energy in units of tx (tx = 1). The bulk of our calculations are for 
ty ~ 1, with tx+y only slightly smaller. This is because broken symmetries other than SC, viz., antiferromagnetism 
(AFM) and CO dominate at weaker frustration^. We first calculate the exact wavefunctions in the lowest total 
spin S = 0 subspace for all p within the periodic 4x4 triangular lattice. In Fig. [^g) we plot the total normalized 
contribution by the covalent VB diagrams with short bonds to the exact wavefunction as a function of p for several 
Hubbard U and V. For moderate to large e-e interactions the maximum in this contribution occurs at p = 0.5, where 
the wavefunction consists predominantly of VB diagrams with short bonds, indicating that VB diagrams with short 
bonds as in Fig.[^f) dominate at p = 0.5. 

In the PEL the occurrence of local singlets is enhanced by correlations and peaks for p ~ 0.5. We anticipate Bose 
condensation of singlet pairs within the PEL state within the mechanism of SC proposed by SchafrotlPil. Without e-p 
coupling in Eq. there is however no static SG and PEC insulating state, which suggests that a complete theory of 
SC will necessarily require explicit inclusions of both e-e and dynamic e-p interactions. As is however well established 
from studies of CDWs and SDWs, the tendency to the dominant instability in models containing both e-e and e-p 
interactions can be determined from studying correlation functions of the electronic Hamiltonian alon^i^. We have 
therefore performed calculations within Eq. [^to determine if the dominance of VB diagrams with short singlet bonds 
at p ~ 0.5 is accompanied by enhanced long-range superconducting pair correlations. The goal is to demonstrate that 
the PEL is a distinct ground-state phase and is a precursor to a correlated superconducting state. 

We define the standard singlet pair-creation operator 
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where giv) determine the pairing symmetry. We have calculated the equal-time pair-pair correlations Pij = (A|Aj), 
for four different periodic anisotropic triangular lattices, 4x4, 6x6, 10x6 and 10x10, with the widest possible carrier 
densities 0 < p < 1, using three different numerical techniques (see Methods). We have chosen lattices that have a 
closed shell configuration at p exactly 0.5, see also Supplementary Information. To facilitate comparison of multiple 
lattices and to mitigate finite-size effects, we calculate the distance dependent pair-pair correlations P(r) (r = \fi — fj\) 
and show here the average long-range pair-pair correlation P = Np X]|r |>2 F’(’")j where Np is the number of terms 
in the surrfi^. 

We have found dx^-y^ and dxy symmetries to dominate over s-wave symmetries in our calculations. Further, for 
each lattice only one of the two d-wave channels is relevant; dx^-y'^ for 4x4 and 10x6, and dxy for 6x6 and 10x10 (see 
Supplementary Information). The origin of this lattice dependence is currently not understood; note, however, that 
the distinction between dx'^-y'i and dxy symmetries is to a large extent semantic in the strongly frustrated regime we 
investigate. Furthermore, it is possible that the actual pairing symmetry is a superposition of dx 2 -y 2 and dxy We have 
not attempted to find this superposition. Rather, for each lattice and p we have calculated the dominant symmetry 
P as a function of U. Plots of P versus U for the different lattices and p are given in the Supplementary Information. 
The complete results, summarized in Figs. and are remarkable: Coulomb interactions enhance superconducting 
pair correlations only in a narrow density range close to p = 0.5. For each lattice P{U)/P{U = 0) > 1 for a single p 
that is either exactly 0.5 or one of two closest carrier fillings with closed shell Fermi level occupancy at U = 0. Pair 
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correlations are suppressed ty [/ at all other p, including the region 0.7 < p < 1 that has been extensively investigated 
in the context of cuprate SCP. In three of four lattices in Fig. ^enhancement of P{U) occurs for p slightly away from 
0.5. The magnitude of pair correlations depend on both the pair binding energy and the kinetic energy to be gained 
from pair delocalization; in finite lattices both quantities depend strongly on the details of the one-electron energy 
spectrum. We show in the Supplementary Information that for each of the four lattices the p at which enhanced P{U) 
occurs can be predicted from the known one-electron levels. Importantly, the deviation from 0.5 of the p at which 
P{U) is enhanced (excluding the 6x6 lattice where this deviation is zero) decreases monotonically with size. 

Having nonzero Vij affects lattice frustration minimally when all three components, 14, Vy and Vx+y are nonzero. 
Pair-correlations for Vx = Vy = Vx+y could be calculated only for the 4x4 lattice, where the behavior of the pair 
correlations is qualitatively similar to Vij = 0, although the magnitude of the enhancement is smaller. We have found 
that when 14 = Vy, Vx+y = 0, dxy pair correlations are enhanced uniquely for p ~ 0.5. Similarly, 14+^ ^ 0, and 
any one of 14, Vy nonzero enhances (suppresses) dx 2 -y 2 {dxy)- Overall, there is a broad parameter region over which 
the pair correlations remain enhanced at the same p where enhancement is found for U ^ Q. (see Supplementary 
Information). 

The numerical ground state results are further conhrmed by DQMC calculations. The sign problem is severe for 
large U, but up to 1/ = 2 the results are reliable even for the largest [3 {(3 = tx/ksT) we have investigated (see 
Methods and Supplementary Information). Fig.shows that with increasing /3 there occurs progressive enhancement 
of P{U) with increasing U, uniquely at p ~ 0.5. As in Fig. U suppresses pair correlations at all other p at large 
(3. The excellent agreement between P{U) obtained from Path Integral Renormalization Group (PIRG) and DQMC 
indicates that while the DQMC calculations could be performed at the smallest T only for U < 2, enhanced pair 
correlations should be expected at even larger U. Fig. [^summarizes the enhancement of pairing as a function of p for 
all lattices, including in the nondominant channels. As seen in Fig. [^ the dominant pairing symmetry is enhanced 
only for p ~ 0.5. Pairing in the non-dominant channels is enhanced weakly for small t/ « 1 for some p, but are 
weakened further as U is increased. 

Our computational results have direct implication for the mechanism of SC in the quasi-2D CTS. Typical supercon¬ 
ducting CTS are the families (BEDT-TTF) 2 X and Z[Pd(dmit) 2 ] 2 - The number of holes (electrons) p per BEDT-TTF 
cation (Pd(dmit )2 anion) is thus 0.5. The same stoichiometry is true for all superconducting (but not merely conduct¬ 
ing) CTS with different organic molec ular components. Strong e-e interaction and frustrated anisotropic triangular 
lattices also are common feature^l^*^. Based on the crystal structures (BEDT-TTE) 2 X are referred to as a, j3, 6, 
K etc. The K -family has been investigated the most intensively, because of the proximity of magnetic phases to SC 
in thes d^^b Sl family, X=Cu[N(CN) 2 ]Cl (k-C 1) is AFM at ambient pressure, X=Cu 2 (CN )3 (k-CN) is a quan¬ 

tum spin liquid (QSL), and X=Cu[N(CN) 2 ]Br (w-Br) and X=Cu(NCS )2 (k-NCS) are superconductors at ambient 
pressurd^^. SC is also observed in k-C 1 and k-CN under pressurd^. Pseudogap behavior has been observed in the 
“metallic” K-Br and k-NCS for T < Tpc ~ 50K from NMRE^M^, STlVpi], precision lattice expansiorP^, and magnetic 
torque measurement^^. 

Structurally, in the K-family the BEDT-TTE molecules are strongly coupled as dimers, which form an anisotropic 
triangular lattice. AEM in k-C 1 and QSL behavior in k-CN are then easily explained if the dimers and not the 
individual molecules are considered as effecti ve site s. This gives an effective p = 1 Hubbard model that will yield 
the AEM (QSL) for weak (strong) frustratior l^'^bsi Precise numerical calculations have however demonstrated the 
absence of SC within the p = 1 Hubbard model for any frustratiorP^^^. Our present work is able to explain both 
the magnetic behavior and SC: in the localized insulating phase the dimerized p = 0.5 and the effective p = 1 model 
give same behavioi!^, but in the pressure-induced delocalized phase larger interdimer hopping leads to breakdown 
of the effective picture and quantum critical transition from the AEM to the PEL, which is superconducting once 
pair coherence is reached. Strong support for this viewpoint is obtained from the recent observation of the PEC in 
K-Hg(SCN) 2 CP. Eurther, our proposed mechanism gives for the first time a cogent explanation of the pseudogap in 
K-Br and k-NCS. The effective p = 1 model fails to explain the T < Tpo NMR behavioil^. Within our theory at 
Tpc there occurs the quantum critical transition to the PEL, while phase coherence and SC are reached only at Tc. 
Diamagnetism observed from Nernst coefficient measurements at T > tJ^ supports this picture of preformed pairs. 
The close proximity of the PEL to the PEC explains the density wave-like behavior in this phase noted by Muller et 
al.l^. In Fig. § enhanced pair correlation at p = 0.5 begins to appear already at /3 = 8; with average \t\ ~ 0.1 eV, we 
see that T pg can be as high as ~ 100 K. 

In other CTS, the semiconducting stateproximate to SC exhibits CO. Pressure-induced SC from CO states is 
seen in a- and 0-(BEDT-TTF) compound^ EtMe 3 P[Pd(dmit) 2 ]J^ and /?-(meso-DMBEDT-TTF) 2 X (X = PFg and 
AsFg)!^. The horizontal stripe CO below the SC transitions in the a- and 0-(BEDT-TTF)p2llI]^ so-called valence 
bond crystal order in EtMe 3 P[Pd(dmit) 2 ] 2 P and the checkerboard CO in /3-(meso-DMBEDT-TTF)2X (X = PFg and 
AsFg)P^all have the same pattern as the PEC in Eig. [^e) (see reference [321 and Supplementary Information). The 
bandwidth-driven quantum critical transition is now directly from the PEC to SC. The strong role of pho nons in SC 
seen experimentalljESl is expected, as it is the co-operative effect between the e-e and e-p interactioni^^ that drives 
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the transitions to the PEC and PEL. 

Two of us have recently pointed out the unusual abundanc^^ of seemingly unrelated correlated-electron materials 
that are superconducting at carrier concentration p ~ 0.5. In all cases the superconductor belongs to a family of 
materials with varying p, and no SC is observed for p substantially different from 0.5 (as is also true for the CTS). 
While the experimental systems should be investigated individually, it is conceivable that the shared features of 
p — 0.5, lattice frustration and strong e-e interaction point to a new paradigm for correlated-electron SC. 
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FIG. 1: Effective e-e attraction in p = 0.5. (a) p = 0.5 dimers with weak (left) and moderate (right) interdimer singlet 
bonding. Sites colored gray, blue and red have charges 0.5, > 0.5 and < 0.5, respectively, (b) and (c) 2kF spin singlet states 
in the p = 0.5 ID chain, for small to intermediate U and V (t/ ~ 4, 1/ ~ 1), and for intermediate to large U and V {U < 10, 
F < 3) respectively. In both cases V < Vc{U^. (d) The PEC in the p = 0.5 zigzag laddei®. (e) The PEC in the anisotropic 
2D triangular lattic^. The CO has pattern .. .1100... in two directions and . . ..1010... in the third direction, where T’ and 
‘0’ denote charge-rich and charge-poor sites. The 2D PEC has the same charge structure as would be obtained from coupled 
zigzag ladders. Double, single and dotted bonds in (b) - (e) denote bonds with decreasing strengths, with the double dotted 
bond weaker than a single bond but stronger than a single dotted bond. Differences in bond strengths result from nonzero e-p 
coupling, (f) Covalent VB diagrams with short bonds in p = 0.5. (g) Exact total normalized weights of covalent VB diagrams 
with short bonds in the ground state wavefunction of different p for the 4x4 lattice, for ty = 1, tx+y = 0.8. 
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FIG. 2: Density dependence of dominant ground state pair-pair correlations. Average long range pair-pair correlation 
P{U) normalized by its uncorrelated value P{U = 0), for the (a) 4x4, (b) 6x6, (c) 10x6 and (d) 10x10 anisotropic triangular 
lattices, for ty = 0.9 and tx+y = 0.8. The 4x4 results are exact; 6x6 and 10x6 results are obtained using the PIRG method; 
and 10x10 by the GPMG method. P{U)/P{U = 0) > 1 for a single p in each case, either for p — 0.5 or for one of two closest 
carrier fillings with closed shell Fermi level occupancies at U = 0. 
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FIG. 3: Temperature dependence of pair-pair correlations, as calculated using DQMC. P{U) for d^y pairing as a 
function of p and inverse temperature P for the 6x6 lattice with G =1, ty =0.9, and tx+y=0-S, calculated using DQMC. Note 
the gradual enhancement of P{U) for p ~ 0.5 with increasing /3, beginning from /? = 8. P{U) is suppressed by U at all other p. 
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FIG. 4: Summary of p-dependence of ground state pair-pair correlation enhancement. Average long range pair-pair 
correlation P{U) normalized by its uncorrelated value P{U = 0), for (a) U = 1 and (b) t/ = 4 (t/ = 2 for the lOx 10 lattice). All 
results with P{U)/P(U = 0) > 1 are included, for both d^ 2 _y 2 and d^y pairing symmetries. The dominant pairing symmetry 
for each lattice is indicated with darker shading. The width of each bar is 1/A^, where N is the number of lattice sites. 



























































